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X-RAY CONSTRAINTS ON THE WARM-HOT INTERGALACTIC MEDIUM 
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ABSTRACT 

Three observational constraints can be placed on a warm-hot intergalactic medium (WHIM) using 
ROSATVSPC pointed and survey data, the emission strength, the energy spectrum, and the fluctuation 
spectrum. The upper limit to the emission strength of the WHIM is 7.5 ± 1.0 keV s~^ cm~^ sr~^ 
keV^^ in the | keV band, an unknown portion of which value may be due to our own Galactic halo. 
The spectral shape of the WHIM emission can be described as thermal emission with logT = 6.42, 
although the true spectrum is more likely to come from a range of temperatures. The values of emission 
strength and spectral shape are in reasonable agreement with hydrodynamical cosniological models. The 
autocorrelation function in the 0.44 keV < E <1.21 keV band range, w{6)^ for the extragalactic soft 
X-ray background (SXRB) which includes both the WHIM and contributions due to point sources, is 
< 0.002 for 10' < 6* < 20' in the | keV band. This value is lower than the Croft et al. (2000) cosniological 
model by a factor of ~ 5, but is still not inconsistent with cosmological models. It is also found that the 
normalization of the extragalactic power law component of the soft X-ray background spectrum must be 
9.5 ± 0.9 keV s'^ cm'^ sr~i keV'^ to be consistent with the ROSAT PAX-Sky Survey. 

Subject headings: cosmology: diffuse radiation — cosmology:observations — X-rays:ISM 



1. INTRODUCTION 

The problem of the missing baryons is well known; mea- 
surements of D in high redshift absorption line systems and 
application of Big Bang Nucleosynthesis models reveal a 
baryon fraction, f^/i^g = 0.0394 ± 0.0029, (Buries & Tytler 
1998) while inventories of local material can account for 
only nh^^-^ = 0.014± 0.0034 (Fukugita et al. 1998). Cos- 
mological hydrodynamical simulations (Cen et al. 1995, 
Gen & Ostriker 1999, Bryan & Norman 1998, Dave et al. 
2000, and Croft et al. 2000 for example) suggest that 
the missing baryons exist as shock-heated low-density fila- 
ments of a "warm-hot" intergalactic medium (WHIM) at 
temperatures of log T = 5 — 7. The visibility of the WHIM 
is compromised by a number of Galactic foregrounds that 
exist in the lower part of this temperature range, the Local 
Hot Bubble (logT = 6.11, Snowden et al. 1998) and the 
Galactic Halo (logT — 6.08 and, perhaps logT = 6.45, 
Kuntz & Snowden 2000a). All of these components emit 
predominately below E=l keV, where Galactic absorption 
can be significant. A component at the high end of the 
WHIM range (e.g., logT = 7) would emit strongly in the 
1.5 keV band. However, Hasinger et al. (1998) has directly 
resolved 70-80% of the fiux in the 0.5-2.0 keV range, leav- 
ing little room for diffuse components. Thus, significant 
diffuse emission at the high end of the WHIM range would 
be difficult to reconcile with observations, while emission 
at the low end of the WHIM range becomes difficult to 
detect. 

Without high resolution, high grasp, non-dispersive X- 
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ray spectrometers, there are three observational tests of 
the WHIM models; the total emission, the energy spec- 
trum, and the fluctuation spectrum. Theoretical calcula- 
tions by Gen & Ostriker (1999) suggest that the WHIM 
produces ~ 7 keV s^^ cm^^ sr^^ keV^^ at 0.7 keV; the 
Croft et al. (2000) prediction for the entire IGM in the 
0.5-2.0 keV band yields the same value for the | keV band 
{ROSATYlAb band). Phillips et al. (2000), however, esti- 
mate a substantially lower X-ray emission, ~ .7 keV 
cm~^ sr~^ keV~^ in the 0.5-2.0 keV band. Dave et al. 
(2000) using a number of models including that of Croft 
et al. (2000) predicts that the baryon fraction as a function 
of temperature will be peaked near logT = 6.45. Finally, 
Croft et al. (2000) predict that the value of the AGE of the 
entire extragalactic SXRB in the 0.5-2.0 keV band should 
be ~ 0.005 for 6 = 10'. 

2. THE EMISSION STRENGTH 

To determine a limit to the strength of the WHIM 
emission, we apply the methods of Kuntz & Snowden 
(2000a) to the north Galactic polar cap, 45° < £ < 270° 
and b > 85°. The left panel of Figure 1 shows broad- 
band data from the i? 05"^ T All- Sky Survey {RASS), for 
the north polar cap as well as the best fit model spec- 
trum, which is shown at a much higher energy resolution 
than measureable with ROSAT. The broad-band data has 
been corrected for band-to-band point source detection 
limit differences. We have assumed that the soft X-ray 
background (SXRB) is composed of the following fore- 
grounds/backgrounds which may be successively removed. 
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Fig. 1. - Left: The Total Soft X-ray Background. The boxes are the RASS broad-band measurements. The soUd lines are our model with 
uncertainty at 100 eV resolution. The dotted line is the contribution due to unresolved Galactic stars after absorption by the Galactic disk. 
The dashed line is the contribution by unresolved extragalactic sources (the extragalctic power law) after absorption. The shaded band is 
the contribution due to the LHB. Right: The Galactic halo emission and the extragalactic diffuse emission. The dotted line is the spatially 
variable soft component that is likely due to the halo. The dashed line is the uniform hard component which may contain contributions from 
both the Galactic halo and extragalactic emission. The bar is the maximum uniform i keV band emission allowed by the measure of Snowden 
et al. (2000). 



Unresolved Galactic Stars: The model of Kuntz & Snow- 
den (2000b) (similar to the model of Guillout et al. (1996)) 
has been used to calculate the spectrum of the Galactic 
stars below the RASS point source detection limit. This 
model reproduces the number of stellar X-ray sources in 
the ID"** to 10^'^ counts ti^^ range found in deep point- 
ings, to about 10%; unfortunately the luminosity function 
is expected to rise for another two orders of magnitude, so 
there may be larger uncertainties for the total unresolved 
flux. 

Unresolved Extragalactic Point Sources: The measured 
spectrum of the unresolved extragalactic point sources has 
been taken from Chen et al. (1997) who find a photon in- 
dex, r = 1.46, from individual RO S AT a.nd ASCA spec- 
tra. Since the RASS has a better non-cosmic background 
subtraction than is possible for any single pointed obser- 
vation, the normalization of this extragalactic power-law 
was adjusted so that all of the known emission compo- 
nents account for all of the band R7 (1.05 kcV < _E < 2.04 
keV) fltix. The required normalization is 9.5 ±0.9 keV s~^ 
cm~^ sr~^ keV~^, the bulk of the uncertainty is due to the 
systematic uncertainty in the R7 band flux, ±3.6 x 10~^ 
counts s^^ arcmin^^. This value is considerably less than 
the recent SAX value of 11.7 ± 0.5 keV s"^ crn^^ sr-1 
keV-i (Vecchi et al. 1999), and is ~ 18% higher than 
the HEAO-1 value of Marshall et al. (1980). It has been 
assumed that the extragalactic power-law does not have a 
break within the i? 05"^ T" energy range as Mushotzky et al. 
(2000) find no break in the 0.3-10.0 keV range for the ag- 
gregate spectrum of a sample of their faintest Chandra 
sources. 

Local Hot Bubble: The strength of the LHB is deter- 
mined from the variation of observed flux with total ab- 
sorbing column in the manner of Snowden et al. (1998) 
and Kuntz & Snowden (2000a). 



The Remainder: The remaining flux contains contribu- 
tions from all components beyond the Galactic absorp- 
tion that are not due to the unresolved extragalactic point 
sources; the Galactic halo and any extragalactic diffuse 
emission. The right panel of Figure 1 shows this remain- 
der after correction for the Galactic absorption. 

As Kuntz & Snowden (2000a) showed that the emission 
from beyond the Galactic absorption cannot be fit by a 
single thermal component, we fit this remainder with two 
Raymond & Smith (1977) model components; the lack of 
independent spectral resolution elements prohibits fitting 
a greater number of components. The fits to the north 
polar cap produce temperatures similar to those found in 
Kuntz & Snowden (2000a) for 6 > 55°, 45° < ^ < 270°, de- 
spite the differences in analysis (subtraction of the compos- 
ite stellar spectrum and the renormalization of the EPL); 
log Ts = 6.08l°:°4 and logT^ = 6A3t'^f,. 

The strength of the soft component varies strongly 
across the sky while the strength of the hard component 
is quite uniform (Kuntz & Snowden 2000a), at least at 
relatively high galactic latitudes (|6| > 30°). The angular 
variation of the soft component suggests a Galactic ori- 
gin, while the uniformity of the hard component suggests 
that it is due either to a Galactic corona in hydrostatic 
equilibrium, or extragalactic emission. There is no way, 
on the basis of i? 05^ T photometry, to separate Galactic 
halo emission from diffuse extragalactic emission. There- 
fore, the upper limit to the diffuse extragalactic emission 
is the strength of the hard component, 7.4 ± 1.0 keV s"""^ 
cm^^ sr^^ keV^^ in the| keV band. Note that the mini- 
mum J keV emission from beyond the Galactic absorption 
is ~ 7 keV s'^ cm'^ sr'^ keV"^ (Snowden et al. 2000), 
consistent with the j keV flux produced by the hard com- 
ponent. 

One might hope to place a limit on the contribution of 
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the Galactic halo by considering the halos of Milky Way 
analogues. MlOl is the closest face-on Milky Way ana- 
logue, though its D25 is about twice that of the Milky Way, 
and its star formation rate as measured by Lfir/D\^ is 
a factor of ^ 5 lower than the Milky Way. The sun is 
~ 0.351)25 from the Galactic center; at an equivalent ra- 
dius in MlOl the upper limit of the | keV surface bright- 
ness due to MlOl (which would include its halo and some 
emission from its disk) is half the strength of our hard com- 
ponent in the same band. However, for NGC 891. which 
has about half the mass of the Milky Way, a similar D25, 
and a similar Lfir/D^^, an observer at ^ 0.35D25 from its 
center would sec a galactic halo that is ^ 9 times brighter 
than our hard component (Bregman & Houk 1997)! The 
relation between the strength of the halo and the stellar 
formation rate is not sufficiently understood to determine 
the fraction of the hard component that might be due to 
our own Galactic halo. 

3. ENERGY SPECTRUM 

The fraction of baryons at each temperature in the Croft 
et al. (2000) model is shown in Dave et al. (2000). The 
baryon fraction at 2; = is peaked at logT ~ 6.45; if the 
baryon fraction were weighted by the emission at (^ach tem- 
perature, the distribution would be more sharply peaked 
and at somewhat higher temperatures. The contribution 
of baryons at higher z is substantially smaller than that of 
baryons at .2 = and is peaked at somewhat lower temper- 
atures (logT ~ 6.05 at z = 1). Given the caveats above, 
the agreement of the temperature of our hard component 
(logT = 6.43) with the "characteristic" temperature pre- 
dicted by the simulations is comforting, but may be com- 
pletely fortuitous if the bulk of the emission is due to the 
Galactic halo. 

4. FLUCTUATION SPECTRUM 

We have determined the autocorrelation function (ACF) 
for the I keV band {ROSAT band R45) from seven mo- 
saics of deep PSPC pointings. The non-cosmic background 
of each PSPC pointing was removed using ESAS (Snow- 
den & Kuntz 1998), all 4ct point sources were removed, 
the image was restricted to the inner 53', and the result- 
ing data were mosaicked (Kuntz & Snowden 1998). The 
remaining non-cosmic background was removed by com- 
paring the mosaics to the RASS . Using the variation of 
transmitted flux with column density, we determined the 
strength of the LHB in the R1L2 band, extrapolated this 
quantity to the R45 band, and removed it from the im- 
age. As the point source detection limit varies across each 
mosaic, the mosaics were corrected using the cxtragalac- 
tic point source log TV — logS* of Hasinger et al. (1998) 
and the Galactic stellar luminosity function of Kuntz & 
Snowden (2000b). Each R45 band mosaic was then cor- 
rected for Galactic absorption. The resultant image con- 
tains contributions from the Galactic halo, the unresolved 
extragalactic point sources, and the extragalactic diffuse 
emission. Once all the corrections were determined, the 
mosaics were reconstructed using only the central 32' ra- 
dius region of each pointing. The smallest scale for which 
one may meaningfully calculate the ACF is the size of the 
largest PSF, 5'. (Although the ROSAT PSF is ~ 0.25' on 
axis, the size of the PSF grows rapidly with distance from 



the optical axis.) 

For each mosaic we calculated the ACF, 

j:{r-{r)){r'-{r'))V^ 

w{e) = 



j:r+j:r') /2 

i V J . 
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where the sum is over all pairs of pixels scpcrated by R 
is the count rate in those pixels, and w is the statistical 
weight (the exposure time) given to those pixel 

We find that there is a significant mosaic- to-mosaic vari- 
ation in the ACF, which was confirmed by variance analy- 
ses [51 /I) of the mosaics. The ACFs of individual mosaics 
were weighted by the mean exposure to calculate the mean 
ACF shown in Figure 2. As the effective point source de- 
tection limit for a mosaic can be much smaller than for in- 
dividual exposures, the ACF contains a contribution from 
residual point sources whose PSF radii range up to 5'; the 
maximum possible contribution due to the residual point 
sources is shown by the dashed line. Thus, the ACF shown 
is the upper limit for scales smaller than 10'. 

The short line labeled "SH" is the ACF determined by 
Soltan & Hasinger (1994) for the ROSATRQl band (0.73 
keV < E < 2.04 keV); over the interval shown it is nearly 
consistent with zero. The dotted lines show the ACFs 
from Chen et al. (1994), corrected from their energy band 
(ROSAT ha.nd R45 + R67) to our band energy (ROSAT 
band R45), assuming that the fluctuations in the R67 band 
are negligible for these scales. The upper line is their mea- 
sured ACF; the lower line is the ACF that they derived by 
fitting the data with the fimction a{6/l')~'^-^ . 

We chose to measure the ACF in the | keV band as 
the hard component discussed in the previous section will 
have a stronger contribution in this band at higher ener- 
gies, but there will not be as great a problem with Galactic 
absorption as there is at lower energies. However, the only 
pubhshed model ACF (Croft et al. 2000) is for the 0.5-2.0 
keV band. Since the Soltan & Hasinger (1994) ACF for 
higher energies is much smaller than that for the | keV 
band, the ACF for the 0.5-2.0 keV band can be estimated 
from our ACF by scaling it by the quantity 

flux in our band 
flux in 0.5-2.0 keV band J ' ^ ' 

which is ~ 0.4. Taking the ACF for the 5' PSF mosaics, 
wc find that the ACF in the 0.5-2.0 kcV band should be 
^ 0.001 for 10' < 61 < 20'. This is about a factor of five 
smaller than the model predictions. 

There may, in fact, be contributions to the ACF due to 
the Galactic halo and the LHB, but these are expected to 
be small. The ACF for emission from the LHB has been 
measured for the 10' < 6* < 20' range and is consistent 
with zero (Kuntz 2000). 

Such a direct comparison of theory and observation is 
not correct. The linear size of the Croft et al. (2000) sim- 
ulation (^ 1°) corresponds to ^ 8 Mpc/i^^ at a distance 
of 500 Mpcft"^ and contains an X-ray emitting cluster of 
galaxies. Since the mean separation of clusters is several 
times larger than 8 Mpcft."'^, the simulation has an atyp- 
ical overdensity (hence atypically large concentration of 
X-ray emission), and thus a higher value of the ACF. 

The observed ACFs show variation among mosaics with 
solid angles of 10 to 30 square degrees; the root of the 
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variance of our set of seven mosaics is about twice the dis- 
played uncertainty. This variation may be consistent with 
cosmic variation. 

5. DISCUSSION 

We have determined the upper hmit to the strength of 
diffuse extragalactic emission in the soft X-ray range, by 
"peehng off" various known foregrounds and backgrounds. 
What is left is some combination of extragalactic diffuse 
emission and Galactic halo emission. We expect our galaxy 
to have an X-ray halo simply due to hydrostatic equilib- 
rium (the virial temperature is ~ 10® K) and because other 
galaxies like our own have them. The strength of the halo 
remains a vexed issue. The combined halo/extragalactic 
emission weighs in at just about the strength expected for 
the extragalactic component alone, so either our galaxy 
has very little halo, or the models over-predict the X-ray 
emission of the WHIM. The lower PhiUips et al. (2000) 
value suggests that more care may be needed in compar- 
ing simulations to measurements. 

It should be noted that this analysis has used a normal- 
ization for the contribution by the unresolved extragalac- 
tic point souces that is at the lower end of the range of 
values commonly used. Independent of the existence of 
extragalactic diffuse emission, or the existence of a Galac- 
tic halo, the RASS will not admit to a higher value for 
the normalization. Increasing the normalization decreases 
the amount of flux available for the extragalactic diffuse 
emission to significantly smaller than that predicted by the 



models. 

The upper limit to the ACF in the | keV band is smaller, 
by about a factor of five, than that predicted by the mod- 
els at a scale of 10', and all the improvements to the mea- 
surement of the ACF (better correction for residual point 
sources and more robust estimates of the contribution to 
the fluctuations by the halo and the LHB) will only fur- 
ther decrease the value of the ACF. If a portion of the 
hard component strength is due to the galactic halo, the 
weakness of the observed ACF signal may be due to the 
reduced strength of the extragalactic emission. However, 
as noted above, one should compare our mean ACF with 
the ACF of a simulation whose size is comparable to the 
cluster correlation length rather than a small region sur- 
rounding a cluster. 

6. CONCLUSION 

The observed emission strength, energy spectrum, and 
fluctuation spectrum of the soft X-ray background orig- 
inating from beyond the Galactic absorption is roughly 
consistent with the predictions by hydrodynamical simu- 
lations of the universe. Even given the caveats concerning 
the Galactic halo, the current observations do not disprove 
the existing models. The observations here present solid 
guideposts for further exploration of the accessible param- 
eter space by future generations of simulations. 

We wish to thank J. P. Ostriker and the anonymous ref- 
eree for their very useful comments. 
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Fig. 2. — The ACF for the j keV band. The ACF has been smoothed by a 3' sUding box. The shaded region is the uncertainty in the 
ACF. The dashed Unc is the maximum possible contribution to the ACF by residual point sources. The dotted lines are the ACF values of 
Chen et al. (1994) corrected from the R47 band to the R45 band. The upper line is their measured ACF; the lower line is their fitted ACF. 
The horizontal line marked "SH" is the Soltan & Hasinger (1994) ACF value for the 1.5 keV band {ROSAT haxid R67). 



